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j>^, Abstract 

^ i' We use covariance and dimensional analysis to find expressions 

for the discontinuity of the potential and normal electric field across 
a flat surface with multipolar charge surface density in vacuum. In 

K*" ' particular we show that 
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0^ . and that 
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'^^- = E E a«,-(5+)^('+"^)(9_)5('— )g,,^(x,y) 
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Mh' Here, E^ is the normal electric field, ^ is the electrostatic potential, 

Qi,m{x,y) is the surface density of the (1^^,171^^) electric multipole and 
k> i d± = d/dx ± id/dy. The prefactors a^^m and 6;^^ in these relations 

^ I are calculated by explicitly evaluating the field of a localized unit 

multipole, both above and below the surface. 
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1 Introduction 

One of the common problems in electrostatics is the determination of electric 
field or potential in the presence of a surface distribution of charges. The 
knowledge of the discontinuity in the electric field across charged surfaces 
is useful in solving boundary value problems in electrostatic. Gauss's law 
allows us to write down a result directly. According to Gauss law the dis- 
continuity of the normal electric field across a surface charge density a is 
47ro" (in CGS units) [|I| . Another important result is the discontinuity of the 
electrostatic potential across a dipolar layer with the dipoles pointing normal 
to the surface. In this case0 5^ = 4:7: P^ where "$ is the potential and P^ is 
the dipole surface density. 

In this paper we study the electric field produced by a flat surface, in vac- 
uum, having general (/*'', m^^) multipole surface density and give expressions 
for the discontinuity in the normal and tangential electric field. We show 
that the discontinuity in the normal electric field is given by 

oo I 



l^^-i V i^ + i V(/ + m)!(/-m)! 

whereas the discontinuity in the potential is and 

oo I I - 

l^^i V i/ + i v/(/ + my.{i - my. 

Here, 

• (li,m{x,y) is the (/*'^,m*'*) multipolar surface charge density over the 
plane. 

• d± = d/dx ± id/dy 

• Ez is the normal electric field 

• \E' the potential. 



From another point of view Eq.(0) may be interpreted as the discontinuity 
in the normal electric field due to an 'effective' surface charge density that is 
given by 
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a^^^ix,y) = J2Y.U^ + (-^y^'^] 
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Similarly, the right hand side Eq.(0) is equivalent to an 'effective' surface 
dipolar layer given by 



oo I 



n^^^(^,y) = EE2-[i-(-iy"^1 ' ^" 



X 



2 Dimensional Analysis & Covariance 

2.1 Definition of the problem. 

Consider a surface at z = charged with (/*'*, m*'*) multiples whose surface 
density is qi^m{x,y). Thus, qi^m{x,y)6x6y is the strength of the {l^^^m*^) 
multipole located at the point (a:, y) on the surface. If ,for example, qi,m{x, y) 
represents a unit multipolar charge at the origin then qi^jn{x,y) = S{x)S{y) 
producing a potential given by 

^_ 47r Yi,m{0,^) 



21 + 1 



4+1 



To fully characterize the surface charge all the multipole density functions 
must be given (recall that / = 0, 1, 2, ... and m = —I, — (/ — 1), ..,(/ — 1), /) 
. For example, for a surface with 'regular' charge density, given by a{x,y), 
the g;,m's are as follows: qofl{x,y) = a{x,y)/^/ATT and qi^mix^y) = for 
higher multiples. These multipole moments give rise to an electric field above 
and below the surface. The purpose of this paper is to find expressions for 
the discontinuity of the normal field and potential in terms of qi^m{x,y). 
Thus, for instance, in the example given above we would expect to find that 
6E.^{x,y) oc qoo{x,y) oc a{x,y) where E.^ is the normal component of the 
electric field. 



2.2 The permissible operators, their transformation 
rules and dimensions. 

In order to find expressions for the discontinuity of the field we use covariance 
and dimensional analysis. We start by stating all the operations that can be 
performed on qi^m{x, y) together with their transformation rule under rotation 
and their dimensions. Using these we show that there is only one possible 
form for the expression, up to an additional multiplicative factor, for the 
discontinuity of the field which is both covariant and dimensionally correct. 

When the system is rotated by an angle a around the z axis we note that 
qi^m transforms according to qi^m -^ e.^^"qi,m (see Appendix A) . The potential 
\1' is a scalar under rotation. Next, any expression for the fields discontinuity 
must depend locally on qi^m{x, y) as the field near the surface (i.e. E(a:, y,z ^ 
0^) ) is determined mostly by the presence of charge distribution at the 
vicinity of the point {x, y,z = 0). 

The first operators that come to mind are multiplication by x and\or y. 
These operators, however, are not permissible since their use would break the 
translational symmetry that is supposed to hold. The only operators that 
are both local and translational invariant are d/dx and d/dy. However, each 
of these operators separately has no definite transformation under rotation 
around the z axis. Operators which do have definite transformation rules 
are easily constructed from these. Consider the operators d± = d/dx ± 
id/dy . When the coordinate system is rotated around the z axis by an 
angle a to form a new coordinate system- {x ,y ) the newly defined operators 
transform according to (see Appendix A) d± -^ e^'^°'d±. Note also that the 
operator 9_9+ is a scalar in this sense and is nothing but the two-dimensional 
Laplacian (which can also be written as d-d+ since d+ and d- commutes). 

We now turn to the transformation rules for the electric field. Obviously, 
the component of the field along the z axis acts as a scalar under the rotation. 
As for the tangential components we again construct new fields defined by 
E± = E^ ± iEy. As E± = —d±'^ and \1' is a scalar the transformation rules 
for E± are the same as those for d±, namely E± — * e^^^E^. 

The following table summarizes the transformation rules under rotation 
for each of the quantities above together with their dimensions: 



Transformation Rule 


Dimensions 


C[l,m -^ e qi^rn 


[Charge] [Length]^^"^ 


d+ -^ e-*"9+ 


[LengthY^ 


Q_ _^ e+iaQ_ 


[LengthY^ 


^ ^^ 


[C harge][LengthY^ 


Ez^Ez 


[Charge] [Length]^"^ 


E+ -^ e-*"E+ 


Charge] Length]^'^ 


E_ -^ e+*"^_ 


Charge] Length]^'^ 


d+d- -^ d+d- 


[Length]'"^ 



2.3 Covariant expression for the discontinuity in the 
normal electric field Ez- 

We start with qi^m{x, y) which transforms hke e*"^". To find the discontinuity 
in Ez{x,y) we operate with d+ a total number of M times and with d- a 
total number of N times to get 

6Ez{x,y)<x{d+)''{d^fqU^,y) (3) 

Applying the transformation rules on both sides of the equation yields the 
first connection between N and M, namely: 



1 p—ictM p+ictN piam 

t 

= -M + N + m 
Dimensional analysis, on the other hand, requires that 

[Charge][Length]~'^ = [Charge][Length]-^-^+^-^ 

t 
-2 = -M -N + 1-2 

Solving for M and A^ gives 

M=\{l + m) N = \{l-m) 

Thus, the general expression for the discontinuity in E^ is 

oo I 

SEz{x,y) = J2Y. aUd+f^'^'^"'\d-P^'^UiA^.y) 

l=Om=-l 



(4) 



(5) 



(6) 



(7) 



where an additional multiplicative factor has been included to make an equal- 
ity out of the proportionality ( Note that the multiplicative factor depends 
on I and m as there is no reason to assume otherwise). The result is further 
summed over all the multiples . 

2.4 Covariant expressions for the discontinuity in the 
potential ^ and the tangential electric field E+ &; 

E_. 

Similar arguments lead to the following result for the discontinuity in \1': 

oo I 

l=Om=-l 

3 Determination of a/^ and 6/^. 

In this section we actually resolve the problem of finding the discontinuity 
of the fields across a multipolar layer. This time however, the exact math- 
ematical expressions are derived together with the required multiplicative 
factor. 

3.1 Fourier Decomposition of Spherical Harmonics Mo- 
ments. 

Let ^i^rn{r,9,ip) denote the potential due to a localized (/,m)-moment with 
unit strength. Thus H], 



21 + 1 r' 



^Ur,g,^) = J; /'ar^ (9) 



The total potential, "^ im{x , y , z) , is then given by a superposition of all the 
contributions from the surface, namely 



where the notation * stands for the two dimensional convolution operation. 



The unit multipole potential may be rewritten in terms of its Fourier 
transform over the {x, y) plane. Since ^im is singular at the origin, the plane 
wave expansion for the upper hemisphere {z > 0) may be different from 
that for the lower hemisphere {z < 0). To make this point clear consider 
for example the case of a dipole at the origin pointing in the z direction. 
Here ^ifi{x,y, Zq) = —^ifi{x,y,—ZQ) so that the expansion of $ for the 
upper hemisphere differs in sign from the expansion for the lower hemisphere. 
Therefore, the plane wave expansion for ^im should generally be written as 

$,^(x, y, z) = ±^jl^l^{k,, ky)e'^'^''+'yy^e-V^^^y\^\dkJky (11) 

where the + sign applies for z > and the — sign for ^ < 0. The dependence 
on z has been included by using V^$;m = and demanding that lim ^im = 

r >oo 

0. A similar expression for qim reads: 

qimix, y) = ^jJQimik., ky)e'^'^^+''yyUkJky (12) 

When Eq.(|lT]) and Eq.(0) are inserted into Eg. (p!0|) , and the identity 

ei(k^^+''yy)dxdy = An^5ik^)Sik^) 

is used, the result becomes 

<i/imix,y,z) = 2^jjqiUk.,ky)^Lik.,ky)e'^''^^+''yy^e-V^^^\^\dk^dky. 

(13) 

In Appendix B we show that ^^{k^, ky), as defined by Eg. ([TTD , is given by: 



^^(k k)- ^ ■,z>0]2n An (-1)- 

''"^ ^' '' ~ I (-1)'+'" ; ^ < / .- V 2/ + 1 V(/ + m)!(/-m)! 

(14) 

{ki + kiY'-^-'^/'{k,+ikyr. 



Substituting Eq.(|T5) for ^f^{K, ky) into Eq.(|TB|) gives 
1 



^im{x,y,z) 



z > I 27r / 47r 



4^2 I (_!)'+- ;z<0 J z-V2/ + l^(/ + m)!(/-m)! 



:-ii 



(15) 



l"'z ~r '''H/ ('''a:: "■" ^"'J/J Qlm[l^xi ^y)^ ^ v ^ \J x y (H^^dl^y 



3.2 The Discontinuity of the potential ^ across the 
Surface. 

The discontinuity of the potential across the surface (denoted by b'^in^ is 
defined as 



b^lra = ^Imix, y,Z = 0+) - ^/^(x, y, Z = Q-). 

Using Eq.(|l^) in Eq. ([T6|) gives 

6%U-,y) = I [1 - (-1)'+™] JJI~ ^-]]J X 

« ' M 21 + 1 ^{l + m)\{l - m)\ 



(16) 



(17) 



47r2 
27r 






[I - (-1)'+™] 



yj yx I "'"yj Himy^x^ '"yj^ u,i\jx'-^'^y 

' 1 \m 



1 / 47r 



2/ + l^(/ + m)!(Z-m)! 



X 



92 92 
+ 



(9x2 ^^2 



(«-m-l)/2 



dx dy 



Qimix,y), 



where in the last equality we have used the fact that a multiplication by 
i{kx ± iky) in {k^, ky)-space is equivalent to the differentiation id± in {x,y)- 
space. Rearranging the differentiation operators in Eq. ([l7D yields 



5<i/,^{x,y) = — \l- (-1)'+"^] -\ 4^ , 

^ ^^ t""^ ^ ^ J ^"^V2/ + l^(/ + m)!(/-m)! 



(-ir 



X 



id. 



i[(/-l)+m](^^_^)i[G-l)-m] 



(!i,m{x,y). 



Comparing this last result with Eq. (^ we first see that the covariant form 
of Eq.([l7|) is identical with what was predicted earlier, and that the prefactor 
bi^rn is given by 



bl^m = 27r [1 - (-1)'+™] 



An 



2l + l^(l + m)\{l-m)\ 



(18) 



For the case of / = l,?7i = for which qifi(x,y) = yS/^vrP^ (a uniform 
double layer with density P^) we find that S'^im = 4:7iPz as required. 



3.3 The Discontinuity of Ez across the Surface. 

The discontinuity of the normal electric field across the surface (denoted by 



6{E, 



iLm 



d 



d 

—'^im{x,y,z = ) 



(19) 



Thus, using Eq.(|T5|) we find that 

97r / Att 

6{E,)U^,y) = -[l + {-iy^^'^] 



, ^ ' = X (20) 

2/ + 1 ^(/ + m)!(/-m)! ^ ^ 



47r2 



-r 



— \l + f-l)'+™l ^ ,/ ^^ 

i""^ J z- V 2/ + 1 ^(/ + m)!(/-m)! 



X 



g2 Q2 

+ 



dx^ dy'^ 



(/-m)/2 



dx dy 



qim{x,y), 



or 



S{E.)U^,y) = 27r [1 + (-1)'+-] ^I^- 



v/(/ + m)!(/-m)! 
(9+)^('+'")(9_)^('-")g,^(x,i/). 



X 



(21) 



Comparing this result with Eq.(^ we see again that the covariant form 
is identical with the one that was predicted earlier, and that the prefactor 
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ai^ = 271 \1 + (-1)'+™] J f"" , ^ (22) 

^ ^ ^ J V 2/ + 1 ^(l + m)\{l-m)\ ^ ' 

Note that ai^m vanishes for odd {l + m) whereas hi^m vanishes for even (/ + m). 
Consequently, for a given (/, m) there can be only discontinuity either in the 
normal field or in the potential but not in both simultaneously. 



4 Summary 

In this paper we derived expressions for the electrostatic field discontinuity 
across a flat surface in vacuum, having general {l^^^m^^) multipole surface 
charge density. We showed that dimensional analysis and the principle of 
covariance alone are suffice to determine the form of the expression. The 
exact expressions for the field discontinuity were then derived mathematically 
and were shown to confirm with our early prediction. These expressions are 

oo I I z ^ 

SE.{x,y) = J2J:^rr[l + i-ir-]J I x (23) 

l^r^-i V 2/ + i^(/ + m)!(/-m)! 

(9+)5('+-)(9_)^('-™)g,,„(x,|/), 
for the discontinuity of the normal electric field, and 

oo I j z _. 

l^^-i V 2/ + i ^(/ + m)\{l - m)! 

for the discontinuity of the potential. Here, qi^m{x,y) is the {p^,m*^) multi- 
polar surface charge density over the plane, d± = d/dx ± id/dy, E^ is the 
normal electric field and \E' is the potential. 

It is interesting to note how these expressions can be extended to the 
electrodynamic case (especially to the discontinuity of the magnetic field). 
To give the reader a sense of the subtleties that arise when considering the 
electrodynamic case we point out that, here, further complication arises as 
there is an additional parameter that must be taken into account, namely /cq, 
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the wave number. The arguments that led us earher to a unique form for the 
discontinuity relation are no longer valid. Instead, we must now incorporate 
the presence of fco into the expressions in a covariant form. This procedure, 
however, is not uniquely determined as is shown by the following argument: 
The operator (i9+)^ must now be replaced by an operator that has similar 
dimension, transforms as (9+)^ and may include k^ as an additional param- 
eter. The general operator that obeys these restrictions is P/v(f^+; ^oe"**^) 
where PN{x,y) is a general two-dimensional polynomial of degree A^ and 
(f is Tarr^iy/x). Similarly, {d-)^ should be replaced by QAr(9_, fcoe^*'^) 
where QN{x,y) is yet another two dimensional polynomial. Covariance and 
dimensional analysis alone cannot yield the coefficients of P/v and Qn- 

Yet another interesting question is how does the discontinuity of the fields 
depends (if any) on the curvature of the surface. Recall that the discontinuity 
in the normal electric field due to a surface charge density, as well as the 
discontinuity in the potential due to a surface dipolar density does not depend 
on the curvature of the surface. One might then carelessly conclude that 
for higher multipolar densities this is the case too. However, the existence 
of the differential operators in Eq . (|23|) , which do not appear at (/,m) = 
(0,0) and (1,0), might suggest that for higher multiples the curvature does 
appear. We conjecture that the curvature does enter in the form of a covariant 
differentiation which replaces the d± operators. 
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A Transformation rule for qi^m{x^y) and d±. 

When the system of coordinate is rotated by an angle a around the z-axis, 
the coordinates transform as 



ry V rtrt /y 

if ^ if' = if — a 

d/dr — > d/dr' = d/dr 

d/dif' = 8/ dip 



(24) 



Thus, while the potential due to a multipolar charge qi^m at the origin is 

47r Yi^^{e,^) 



^i,m{r,0,ip) = qi. 



'21 + 1 r'+i 



as expressed with the old coordinate system, the potential in the new coor- 
dinate system is 



^'Lmir^e'.ip') =qi. 



47r Yi,^{e\^' 



'•'"2/ + 1 r"+i 
Since they must give the same value we have 

qi,mYi,m{e, if) = q'i„Xi^m{0', ifi') = q'i,mYi,m{0, if - a) 
Using the properties of ^z,m() which depends on ip only through e*"**^ we have 

which shows that g[ ^ must be given by 

r,' — r, „ima 
Hl,m, — qi,m(i 

We now turn to find the transformation rule for d±. Rewriting d± in terms 
of polar coordinate system yields 



OX oy 
When these are used with Eq.(^4l) one finds that 



a 




i a 1 




+ 




Or 




r dip_ 



(25) 



dL = e^''' e ' 



-'± 



did 
dr r dip 



e^'°54 
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B Derivation of ^^^{kx, ky). 

In this section we show that ^f^{kx.,ky) as given by Eq.(]T5[) is indeed the 
Fourier transform of <l>/m(a:, ?/, 2; = 0). The Fourier transform of $;„, over the 
x-y plane may be calculated straighforwardly by using the following definite 
integral PI 



00 





,itJ--'2p-ay 



xia' + x^)-^^p-^\ [^^J Uxy)dx = %^ 

{Re{a) > 0, y > 0, Re{v) > -1, Re{fi) > |). 

However we prefer to show it by explicitly evaluating the inverse Fourier 
transform. 

Substitution of Eq. (|^) inside Eq. (|ll]) gives: 



1 ;z>o 1 1 r^^ i-ir 

X 



(-1)'+™ ;z<0 j 27r2™ V 2/ + 1 ^(/ + m)\{l - m)\ 
{kl + klf-'^-^^/^ik, + tkyre'^''-''+^yy^e-Vi^^^\'\dk^dky. 

By using polar coordinates (defined by kr = ^ykl + k"^ and ipk = Tg~^{ky/kx)) 
Eq. (p^) becomes: 



1 ;^>o 1 / 477 (-ir .271 

(-iy+- ;2;<0 i 27rz™V2/ + l^(/ + m)!(/-m)! ^ ^ 

27r 



The integral over the angular part is easily evaluated by using 



27r 

Jm{x) = ^J e*"^°^'^-^'"^rf¥^. (28) 
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This gives: 



(-iy+™ ■,z<0 j 27r^-V2/ + V(/ + m)!(/-m)! ^ ^ 





Since 11 






f x''~^J,{px)e-''''dx = (a" + /32)-/^/2r(/x + i.)p-_^,(a/v/^^^T^) (30) 



a,/5 > ■,Re{n + u) > 0, 



we get 



1 ;^>0 1 1 / 47r (-ir .^.. 

:-!)'+- ■,z<0 j 27rz-V2/ + l^(/ + m)!(Z-m)! ^ ^ 

gim<^27ri"r-('+i)r(/ + m + l)Pf ™(|cos^|) 



By using 



Prix) = i-^y^jY^pn^) (32) 



and the definition of T{x) by which r(/ + m + l) = (/ + m)! we find that 



1 \Z>0 ] I I 4:71 (-1 



' ^ ' ' ' X (33) 



,im<^27ri'"r-('+i)(/ + m)!(-l)'" ^~"^ ; P"(|cosg|) 

(/ + my. 



rriji 



Using the properties of P]^{x) under inversion it may be easily verified that 
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(_/y+„^ :^<o }pr{\cose\) = pr{cose), m 



and since p 



YUO,^) = W^fq:^^r(cos^)e^--, (35) 



we finally find, as required, that 



Air Yirni e,<f) 
21 + 1 r' 



I= J:/u,^ ='^^rn{rl (36) 
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